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£SJ , Abstract 

The object of this paper is to characterize the third order moments (cumulants) and 
bispectra of a homogeneous isotropic field defined on a plane. We establish a one to one 
correspondence between the third order cumulants and the bispectra of such a process in 
terms of Bcsscl functions. 

Dedicated to Lajos Tamdssy on the occasion of his 90th birthday. 

1 Introduction 

In many real applications associated with random fields, the assumption of Gaussianity is un- 
realistic. For example, consider the data of cosmic microwave background (CMB) anisotropies 
provided by NASA which is a classic case of non-Gaussianity. |Mar04] . CMB data are given on 

^sO • a surface of a sphere. The assumption of non-Gaussianity is needed in the investigation of the 

primordial bispectrum given on the whole space, see [VWHK00] . [YKW07] . [AC12J . In time 

^ \ series analysis the non-Gaussianity has been well studied [SRG80] . [SRG84J, [Hin82j, |TM98j . 



|Ter99j . It is known that for a Gaussian time series the bispectrum and all higher order spec- 



tra are zero, and equally well known is the fact that for a non Gaussian process defined on 
a real line the higher order cumulant spectra and higher order cumulants uniquely determine 
each other [BriOlj . similar questions has been considered recently for fields on spheres |MP10j . 
[Terl3j . No such results are available and well known for a homogeneous isotropic field defined 
/\i . on a plane. The Gaussian homogeneous and isotropic fields are well studied see |Yag87| and 

references therein. 

In this paper we first consider a homogeneous field and define the homogeneous and isotropic 
fields for general not necessarily Gaussian fields. The second order properties are given which 
serves as an example of the method we use for obtaining an expression for the third order 
spectra . The main object of this paper is to describe the unique relation between the third 
order covariances (bicovariances) and the bispectrum of a homogeneous and isotropic field on 
the plane. 

1.1 Non-Gaussian homogeneous and isotropic field on the plane 

We consider a homogeneous real valued stochastic field X (x) on M 2 with X (x) = 0. Let us 
suppose that X (x) is continuous (in mean square sense), its spectral representation is 

X(x)= f e^Z (du) , lv,x£R 2 , 



with a finite spectral measure 

E\Z(du)\ 2 = F (du). 

By homogeneity we mean (in strict sense ) the distribution of X (x) is translation invariant and 
all necessary moments exist, see Yag87| for details. Rewrite X (x) in term of polar coordinates 



/"GO / > 27T 

X(r,(p)= / e iprcoa ^-^Z(pdpdrt) 

Jo Jo 

oo 



' — ' /»oo 

= £ ^ '/ Ji(pr)Z e (pdp) (1) 

where Je denotes the Bessel function of the first kind, x = (r,(f), u= (p, rj) are polar coordinates, 



r = \x\ = \Jx\ + x?>, p = |w|, x ■ u = rp cos (ip — rj) and in the above we used the following 
notation 

/•27T 

Z t (pdp) = / i e e~ iiv Z {pdpdri) . 
Jo 

The representation ([T]) will be an orthogonal representation if we assume that Fq (duf) is isotropic 

(invariant under rotations) Fq (cLj) = E\Z(diu)\ = E\Z(pdpdrj)\ = F (pdp) dr]. Under this 

assumption we consider the Jacobi- Anger expansion 

oo 

iprcos(^) = ^ i £ J £ (pr)e ie ^- v) , (2) 



e 1 



and substitute this into the spectral representation of X (r, <p). The stochastic spectral measures 
Zg (pdp) defined above has the following property 



r2ir 

Z_ £ (pdp) = / r e e~ ie71 Z (pdpdrj) 
Jo 

r 2TT 

i l e Uri Z (pdpdri) 



(-iyz e (pdp), 



'-Z (dui) = e l -—Z (du) , and this follows from the fact that the field X (x) is real 
valued. Moreover it is orthogonal 

Cov(Z £l (pid Pl ),Z h (p 2 dp 2 )) = I ^ i^-^e-^-^Fo (d V pdp) 

Jo 

= 5 h -i 2 2TiF(pdp) , 

where ^i x -i 2 is the Kronecker delta. The well known property of Bessel functions Ji = (—1) J-e 
implies that in the representation (pQ) the £ term and the — £ term are conjugate to each other. 
Note that the spectral measure of the stochastic measure Zp (pdp) does not depend on £. We shall 
assume in particular cases that F (pdp) is absolutely continuos, i.e. F (pdp) = a 2 \A(p)\ pdp, 
here a 2 \A(p)\ is usually known as the second order spectrum. In other words the representation 
(PQ) has the particular form 



00 /"OO 

X (r, ip) == Y, ^ \ Je (Pr) ^ (p) W e (pdp) 

i=-oo J ° 



such that Wg (pdp) is a series of white noise measures with constant spectrum, hence 

Cov (W tl (pidpi) , W i2 (pzdpi)) = S ei ^ i2 a 2 pdp. 

The general theorem of Yadrenko ([Yad83j Theorem 1. pp.5) on the spectral representation of a 
homogeneous and isotropic field X (r, ip) gives the representation ([T|) with real valued stochastic 
spectral measures Ze (pdp) constructed directly form the field X (r, ip) itself. 
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2 Isotropy on the plane 

We consider rotations about the origin of the coordinate system. Under a rotation (passive) 
g £ SO (2), by which we mean a rotation when vectors remain fixed, but the point it defines 
is given by a new set of coordinates. A rotation g is characterized by an angle 7 and by the 

rotation matrix 

cos 7 — sin 7 

sin 7 cos 7 

If x £ M. 2 is given in polar coordinates x = (r, tp), then gx = (r, ip — 7), and as usual the operator 
A (g) acts on functions / (r, p), such that A (g) f (r, ip) = f (g^ 1 (r, p)) = / (r, p + 7). 

The isotropy usually is defined through the invariance of the covariance structure. This is 
satisfactory for Gaussian cases but for non-Gaussian fields we need invariance of higher order 
cumulants as well. We need a stronger definition to achieve a similar invariance to be able to 
define third order spectra, which we will propose below. 

Definition 1 A homogeneous stochastic field X (x) is strictly isotropic if all finite dimensional 
distributions of X (x) are invariant under rotations. 

If the homogeneous field X (x) is Gaussian, then the isotropy of the spectral measure Fq (du), 
i.e. in polar coordinates Fq (du) = F (pdp) dnj (2vr), implies 

Cov (A (g) X (xj ,A(g)X (x 2 )) = Cov (X (xj , X (x 2 )) , 

for each x 1; x 2 and for every g G SO (2). That is the distribution of a Gaussian isotropic field 
is invariant under rotations. The definition given above is a generalization of this property 
which will be used later. If the field X (x) is non-Gaussian and if we assume the existence of all 
moments then isotropy follows and followed by that all higher order moments are also invariant 
under rotation. Let us consider a homogeneous and isotropic stochastic field X (x) = X (r, p), 
(r > 0, ip € [0, 2ir)) on the plane 



X 



X(r,<p)= V e u r I J e (pr) Z e (pdp) 



00 



where Zi is an array of measures, orthogonal to each other 

Cov (Z £l (pidpi) , Z i2 (p2dp 2 )) = 6 ei _ i2 2nF (pdp) . 

In this representation e v plays a role of spherical harmonics of degree £ with complex values 
on the plane. In this way an isotropic random field X (x) can be decomposed into a countable 
number of mutually uncorrelated spectral measures with a one dimensional parameter. 

The rotation g takes effect on the 'spherical harmonics' e LLp (m = ±1), as A.(g) e LLp = 
e iem(ip+j) _ e iim-y e i£mip ^ s [ nce fa e 'spherical harmonics' e timip can be considered as a function of 
p. The isotropy of 

X(r,p)= Y, <& Je(pr)Z e (pdp), 

implies that the distribution of X (r, if) does not change under rotations g G SO (2). Consider 

00 POO 

A (g) X (r, p) = Y, z UW+l) / h {(»■) ^ (pdp) 
= V e"M J e (pr)e^Z e (pdp) 

00 POO 

= Y ^ '/ Je( P r)Z e (pdp), 



hence the distribution of Zg (pdp) and e % ^Zi (pdp) should be the same. Now it is evident that 
for a Gaussian random field X (r, <p) , the necessary and sufficient condition of isotropy is the 
Zi {pdp) be independent measures. Indeed under isotropy assumption we have 

Cum (Z h ( Pl d Pl ) , Z l2 (p 2 dp 2 )) = e^ tl+ta ^ Cum (Z h ( Pl d Pl ) , Z< a (p 2 dp 2 )) , 

for each 7, hence either t\ + 1 2 = 0, or Cum (Z^ (pidpi) , Zn 2 (p 2 dp 2 )) = 0, therefore 

Cov (Z £l (pidpi) , Zi 2 (p 2 dp 2 )) = 0, 

unless t\ = £ 2 . 

In general, we have under isotropy assumption 

Cum (Z tl {pxdpx) ,...,Z ip (p P dp p )) = e *(*+*»-+«ph Cum (Z tl (pxdp x ) ,...,Z ip ( Pp d Pp )) , 

that is if £i+£ 2 ■ ■ -+£ p = 0, or Cum (Z^ (p\dp\) , Zi 2 (p 2 dp 2 ) , . . . , Z^ p (p p dp p )) = should be ful- 
filled. In turn, if this assumption is satisfied then the cumulants Cum (Z^ (p\dp\) , . . . , Ze (p p dppj) 
are invariant under rotations and if in addition the distributions of X (r, <p) are determined by 
the moments then the field is isotropic. 

2.1 Spectrum 

Briefly, we consider the second order structure of a homogeneous and isotropic stochastic field 
which is well studied for Gaussian cases. For notational convenience let us denote the stochastic 
process 



POO 

zti (r) = Je (pr) Z e (pdp) 

Jo 



and consider the covariance 



(£1=00 £2=00 

£1=— 00 £2=— 00 

/■OO °° 

= 2tt / V e^-riJt ( P n) Je (pr 2 ) F (pdp) 
Jo £=-00 

/•oo 

= 2tt / Jo (pr) F (pdp) , 
Jo 

where r = \x — y .In arriving at the above formula we used the addition formula of Bessel 
functions, see [EMUT54] Tom2, Ch7, 7.6.2.(6), [Yad83j . Now one may derive the same result 
using the properties of homogeneity and isotropy. We have in general 



and in a special case 



x /'CXI 

X(r,tp) = V e^ / J e (pr)Z e (pdp) 
£=-00 Jo 

x /*oo 

X (rN) = V i l J e (pr) Zi (pdp) , 

._ „ Jo 



(3) 



£=-00 

X (0) = / Z (du) (4) 

00 

Z (pdp) , 




where the 'North pole' N = (0, 1). Let , r = \x — y \, C 2 (r) = Cov [X (x) ,X (y)), and use the 
invariance under translation and rotation to obtain 

Cov (X (x) ,X(y))= Cov (X(x-y),X (0)) 
= Cov(X(rN),X(0)) 



/>oo 

2vr / J (pr) F (pdp) . 
Jo 



The above shows one to one correspondence between the second order covariance and its spectral 
density function, see [Yad83| . In particular for absolutely continuos spectral measure F (pdp) = 
a 2 \A (p) | pdp we have 

roc 

C 2 (r) = 27r J (pr)a 2 \A(p)\ 2 pdp, 
Jo 

in turn 

1 f °° 
a 2 \A(p)\ 2 = —J J (pr)C 2 (r)rdr, 

in case that the integral exists, see |Bri74| . [Yag87|, for instance. The above property of Hankel 
transform used above is based on the following property of Bessel functions 

no 

Jo (pr) Jo (kt) rdr = , 



where 5 (p — n) denotes the Dirac 'function', more precisely 5 (p — n) is a distribution, see 
[AWOlj Sect 11. 



3 Bispectrum 

The third order structure of a homogeneous and isotropic stochastic field X (x) is described 
by either the third order cumulants (bicovariances) in spatial domain or the bispectrum in 
frequency domain. Using the spectral representations, we obtain the third order cumulants 
(central moments) and it is given by 

Cum(Xfe),X(x 2 ),Xfe)) = [[[ e ! (^^)s3(y 1 ,y 2 ,y 3 )4(sk) f[dw k , 
under isotropy for each g £ SO (2) 



fc=i 

! xR 2 xl 2 



Cum(X(gx 1 ),X(gx 2 ),X(gx 3 ))= Iff e < s ?&0 S 3 (g^, gu 2 , gu^) 6 (E?u; fe ) ]l^ fc 



3 

n 
fc=i 

R 2 xK 2 xR 2 

= Cum(X(x 1 ),X(x 2 ),X(x 3 )) 

hence S3 (guii, gto 2 , gu 3 ) = S3 (c^,^,^) = S3 (pi,p 2 ,P3) ■ Now we apply the invariance of the 
bicovariance 

Cum (X (x x ) , X (x 2 ) , X (sg)) = Cum (X (0) , X (x 2 - x t ) , X (x 3 - x x )) 

= Cum(X(0),X(\x 2 -x 1 \N),X(g(x 3 -x 1 ))) 



where g denotes the rotation carrying x 2 ~ 2±i m ^° the North pole. We use the particular 
representations dU and Q and get 

oo 

Cum (X (0) , X (r 2 N) , X (x 3 )) = £ Jjj &e ie ^J h (p 2 r 2 ) J h (p 3 r 3 ) 

x Cum (Z (pidpi) , Z i2 (p2dp 2 ) , Z ia (p 3 dp 3 )) 

oo 

= E fff ^JdP2T2)J-dP^) 

£=—oo q 

x Cum (Z Q (pidpi) , Z e (p 2 dp 2 ) , Z_ £ {p 3 dp 3 )) 

where <p = it/ 2 — <p 3 , is the angle between N and x 3 . The third order cumulant of the stochastic 
spectral measure Z (do/) of the homogeneous field X (x) is given by 

Cum (Z {duty) , Z (dw 2 ) , Z (dw 3 )) = 5 (£fu; fe ) S 3 (w 1; cj 2 , w 3 ) du 1 dw 2 duj 3 

3 
fe=l 

where cD fc = co k / \io k \. The stochastic spectral measures Zi (pdp) are connected to Z (da/) by 

/>2tt 

Z e (pdp) = / t / e -iftj Z (drfpdp) , 

JO 

therefore 

Cum (Z (pidpi) , Z e (p 2 dp 2 ) , Z_ e {p 3 dp 3 )) = S 3 (pi,p 2 ,p 3 ) (5) 

2tt 

x JJje- ie ^-^ 8(Y>\p k u k ) Up k dp k d Vk . 

o 

In order to understand the applicability of the Dirac 'function' in polar coordinates we express 
it by an integral through the Jacobi- Anger expansion ([2|) and obtain 

5(x!p k Q k )=^ 2 -JJe<±^)d\ (6) 



1 />0O fl-K 3 OO 

ml / II E i ek Ji k (Pk\)e i£ ^-0 XdXdC 
!vr) Jo Jo t=ie k t^oo 



1-2-iv 3 °° 

, / / 

(2tt)< 

Now substitute. © into[5j The 'spherical harmonics' are orthogonal 



2w 

I" [ f [ e-tto-** Ue li ^-^d m dC = 5 il 5 i2+e 5e. i . i (27T) 4 
Jo JJJ fc=l 



hence obtain from the equations 

2tt 

P /'OO 

e -tf(i»-w0 5 (E?p fc w fe ) d Vl dri2d V3 = (2vr) 2 / J (piA) J_ € (p 2 A) J, (p 3 A) AdA. 

7 JO 



This integral can be evaluated, if \p2 — p 3 \ < p\ < p 2 + p 3 , and let us denote R 
{pI + p\~ Pi) I ( 2 P2P3), p\ = P 2 2 +P 2 3 - 2P2P3 cos (r?), then 

t ( \\ t l \\ t l \\\a\ cos(larccos(#)) 
Jo (piA) J e (p 2 \) J £ (p 3 X) XdX = , 

o 7rp 2 p 3 V± - R z 

6 



see |PBM86j Tom. II, 2.12.41.16, hence 



2tt 




-i/fni-Ttal <- /^f - \ , , / lX / cos (£arccos (i?)) 
e *(* *) ^ ( s 3 pfc ^) dr]ldmdm = 4vr (_!)< v ^ 

P2P3V 1 - 1^ 

Since J^ = (—1) J-t we have 

Cum (Z (pidpi) , Z £ (p2dp 2 ) , Z_ t (p 3 dp 3 )) 

a 1 -i\l n a n cos(£arccos(i?)) 3 

= 4vr (-l)*<5(pA) i 7 =^=^S 3 (pi,P2,p 3 ) EI MP* 

P2Pzvl-R z fe=i 

where 5 (pA) = 6 (p\ + p| — 2p 2 p 3 cosr] — pf) and the wave numbers pi, p 2 , and ,93 should 
satisfy the triangle relation. Now we arrive at the following expression 



Cum (X (0) , X (r 2 N) , X (x 3 )) = 4^ J^ / / / e"^ J, (p 2 r 2 ) J_< (p 3 r 3 



1 -,\l 1 / a \ cos(^arccos(i?)) 3 

x (-1)* 5 (pA) — ^ — — y s 3 (pi,p2,ps) n pa 

P2P3V 1 - -R^ fc=i 



00 

X „ „ r p 2 +P3 



47r E // / e- i ^J e (p 2 r 2 )Ji(p3r 3 ) 



-OO " - -'IP2-P3| 



cos(£arccos (i?)) 3 

-03{Pl,P2,P3) [[Pkdpk 



P2P3 



VY - W -ov^^.o,^ 



£= — 00 Q 

COS (c?7) . . 

x /, o z ^' 5 ' 3 (pi, P2, /03) sin ryd?7 Jl Pfcdp fe 

y / l-COS / (77J fc=2 

00 °° 
= 4vr E // fe'^JeiPz^Jiipzrs) 

£= — 00 Q 

3 

x cos (£rj) S 3 (px,p2,p 3 ) dr] J] p k dp k , 

k=2 



where p\ = p\ + p\ — 2p 2 p 3 cos (77) and 

Pi dpi 

drj 

The function 



p 2 p 3 smr). 



T{(p,p 2 ,p 3 \r],r 2 ,r 3 ) = ^J e ^J f (p 2 r 2 ) J? (p 3 r 3 ) cos (&/) 

£=-oo 
00 

= ^J cos (£ip) Jg (p 2 r 2 ) Je {p 3 r 3 ) cos (£rj) 

e=-oo 

00 

= d (p 2 r 2 ) J (p 3 r 3 ) + 2 ^ cos (^p) J^ (p 2 r 2 ) J e (p 3 r 3 ) cos (#7) , (7) 



gives the transformation of the bispectrum S3 (pi, P2, P3) into the bicovariance C 3 (ri, r 2 , r 3 ) since 
both angles if and ij together with two sides define the third sides p\ and t\ of the triangles, 
given by wave numbers (pi, p2, P3) and distances (ri,r 2 , r 3 ). Notice that the transformation T 
can be simplified 

T(ip,p2,P3\r],r 2 ,r 3 ) = - I J (p 2 r 2 ) Jo (^3^3) + 2 ^ J? (p 2 r 2 ) J^ (p?,r 3 ) cos (£ (rj + <p)) J 

+ 2 ( J ° (P^l) J (P3 r 3) +2^2,J e (P2T2) Jt (P3^3) COS (£ (f) - if)) J 

= 2 ( J o( w +) + Jo(w-)), 



where w+ = \J (p 2 r 2 f + (P3^) 2 - 2p 2 r 2 p 3 r 3 cos ((77 + <p)), 

«•- \Kp2r2) 2 + (P3r 3 ) 2 - 2p 2 r 2 p 3 r 3 cos ((7/ - if)), see |EM()T54j Tom.2, Ch7.7.15. for 



Jo (z) Jo (Z) + 2 JT J^ (z) J e (Z) cos (A?) = J 






(./' 



V-2 2 + Z 2 - 2zZ cos (■&) 



^1 +^ 2 + <^ 3 = 0) hence the wave numbers pi, P2, and p 3 satisfy the triangle relation (they 
should be able to form a triangle). 

We have now established a relationship between the bicovariance and the bispectrum. 

00 

Cum(X(0),X(r 2 N),X(x 3 )) = 2n If I (J (w+) + J (w-)) S 3 (pi, p 2 , p 3 ) drift p k dp k . 

J J JO k=2 



We show that the bispectrum can be also expressed by the bicovariance in turn. The bicovariance 
Cum (X (0) , X (x 2 ) , X (x 3 )) depends on the length of vectors x 2 , x 3 and the angle p> between 
them, this way a triangle with vertices 0, x 2 , x 3 is formed with length of the third side n, such 
that r\ = r\ + t\ — 2r 2 r 3 cos (p). According to this definition of r\, we introduce the notation 

C(n,r 2 ,r 3 ) = Cwn(X(0),X(x 2 ),X(x 3 )). 

We show that 

00 

Ss(pi,p2,P3) = — T3 // / (JoO+) + Jo(w-))C(r 1 ,r 2 ,r 3 )dipY[r k dr k , 

(2-K) J J Jo k=2 



where the earlier notations are used. 
Consider the integral 



I(Pi,P2,Pz\ki,K2,k 3 )= // / I J (P2T2) Jo (fi3r 3 ) + 2 ^ cos (tip) Jt (p 2 r 2 ) Jg (p 3 r 3 ) cos (trj) 

x Jo (K 2 r 2 ) J (n 3 r 3 ) + 2 ^ cos (£<p) J e (k 2 t 2 ) J e (n 3 r 3 ) cos (M) dp ]J r k dr k , 



3 

E 

k=2 

notice first that cos (£ip) is an orthogonal system on [0,7r], i.e. 



/ cos (tup) cos (£ 2 p) dp = 5 £l= i 2 J J . f * 



hence the integral by cp gives 
I(p 1 ,p 2 ,P3\Ki,K 2 ,K 3 ) = ir I J (p 2 r 2 ) Jo («2r 2 ) r 2 dr 2 I Jo (p3r 3 ) J (K 3 r 3 ) r 3 dr 3 



00 00 

00 ^ ^ 

2tt Y^ cos (£tj) cos (£•&) / J e (p 2 r 2 ) J e (K 2 r 2 ) r 2 dr 2 / J e (p 3 r 3 ) J e (n 3 r 3 ) r 3 dr 3 . 





The integral of Bessel functions provides the Dirac function 

J 1 (pr) Jt (nr) rdr = . 



see |AW01| Sect 11. therefore 



00 



I(pi,p 2 ,p 3 \Ki,K 2 ,K 3 ) = 7T 2 5(p 2 ~ K 2 )5(p 3 - K 3 ) - + - V] COS (llj) COS (M) 

\ 7T 7T ' — ' / 

V 1=1 / 

= 7T 2 (5 (p 2 - «2) ^ (P3 - « 3 ) ^ (77 - 1?) , 

since the - and Ty/ir/2 cos (£r]) forms an orthonormal system, see |AW01j Sect. 11. 

Theorem 2 Ze£ us assume that both integral exists below then 

00 

r r i-k 3 

C(n,r 2 ,r 3 ) = 2tt / (J (w+) + J (w_)) 53 (j>i,f>2,P3) df] \[ p k dp k , 

J J JO k=2 



00 

S 3 (pi,P2,P3) = - — 7b- / (Jo(w + ) + Jo(w-))C(r 1 ,r 2 ,r 3 )dipl\r k dr k . 

l2lT) J J Jo k=2 



Note that in some cases it is more convenient to use the transformation ([7]) to obtain 

00 

f f f n 3 

C{n,r 2 ,r 3 ) = 4vr // / T {if, p 2 , p 3 \r],r 2j r 3 ) S 3 {p u p 2 , p 3 ) dr]\[p k dp k , (8) 

J J Jo k=2 



00 

S 3 (pi,p 2 ,p 3 ) = 4^3/// ^(P>P2,P3h,r2,r 3 )C(ri,r2,r 3 )#nrkdrfc. (9) 



A very useful model with finite many parameters on the plane is connected to the Laplacian 
operator. 

Example 3 The spatial white noise dW (r, ip), on the plane is given as a generalized field by 
the series representation 



x />oo 

dW (r, if) = V e^ / J, (pr) W* (pdp) , 



fc-00 



acian 



where Wi(pdp) with E\Wt(pdp)\ = a 2 pdp, see \Yag81^ , }Yad83^ . We consider the Lapl 
field on the plane by the equation 

(V 2 -c 2 )X{r,ip)=dW(r,ip), (10) 



where 



„ 2 1 d d 2 Id 2 



r dr dr 2 r 2 dtp 2 ' 
is the Laplacian operator in terms of spherical coordinates. Now we have 



/"OO PC 

V 2 e^ / J e (pr) Z e (pdp) = e**<> / 
Jo Jo 



dr 2 r dr r 2 



Jt (pr) Zi (pdp) 



POO 

-e U * / p 2 J e (pr) Z t (pdp) , 
Jo 



hence 

POO poo 

(V 2 - c 2 ) j* / J e (pr) Z e (pdp) = -e^ / (p 2 + c 2 ) J e (pr) Z e (pdp) 
Jo Jo 

Let us compare the two sides of the equation [W\) , and obtain 

- (p 2 + c 2 ) Z e (pdp) = W e (pdp) , 
(p 2 + c 2 ) 2 2-kF (pdp) = a 2 pdp, 

F(pdp)= a2pdp 2 . 
^ " 2vr (p 2 + c 2 ) 2 

We have got the covariance 

POD 

Cov (X (x) ,X(y))=27r Jo (pr ) F (pdp) 

Jo 

a 2 pdp 

J (P r ) , 2 , 2 .2 
(p l + C z ) 

,rK-i (cr) 
~2c ' 



a 



This covariance is of Matern Class, see \Whi5$ (K-\ = K\), in terms of modified Bessel 
(Hankel) function. The Bispectrum of the Laplacian model follows 



3 a 2 



S3(pi,p2,ps) = n . , ■■>■ 

k=1 p% + C 2 

where p\ = p\ + p\ — 2p 2 p-$ cos (rj). We apply the transformation ([?[) 

oo 

T(pi,P2,P3\k 1 ,k 2 ,k 3 ) = J (p 2 r 2 ) Jo (P3r 3 ) + 2^ cos {tip) J e (/9 2 r 2 ) J e (p 3 r 3 ) cos (tri) , 

e=i 

for establishing the connection to the cumulants and simplifying the bicovariance. If (a/b) < 1, 
we have 



cos (trj) 



-dr] 



TT 



Vb 2 -a 2 -b 



o b + a cos (n) y/b 2 - a 2 



see IGROOj/ 3.613.1. Put 



a = -2p 2 p 3 , b = pl + pl + c 2 



b 2 -a 2 = (( P2 - P3? + c 2 )(( P2 + P3? + c 2 ) 



hence 



cos (trj) 



;drj 



■k / \/b 2 — a 2 — b 



Pl +P 2 ~ 2 P2P3 COS (ri) + C 2 ~'' yjb 2 - O 2 
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In this way we arrive to the Fourier expansion of the bicovariances 

00 

C(ri,r 2 ,r 3 ) = /o + 2^/icos(ty), 

with coefficients 



/ \ £ 

-ir I v b 2 — a 2 — b\ _ , . _ . . .5. 



/j = 4?r / / 2 1 — - — i j^ (p 2 r 2 ) j£ (p 3 r 3 ) n ^rfpfc- 

^ ' 
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